Spin-lattice coupling, frustration and magnetic order in multiferroic RMnOs 



X. Fabreges^, S. Petit^, I. Mirebeau^, S. Pailhes^, L. Pinsard^, A. Forget"^, M. T. Fernandez-Diaz^ and F. Porcher^ 
^ Laboratoire Leon Brillouin, CEA-CNRS, CE-Saclay, 91191 Gtf-sur-Yvette, France 
^Lahoratoire de Physico-Chimie de I'Etat Solide,ICMMO, Universite Paris-Sud, 91405 Orsay, France 
'^Service de Physique de I'Etat Condense, CEA-CNRS, CE-Saclay, 91191 Gif-Sur-Yvette, France and 
'^Institut Laiie Langevin, 6 rue Jules Horowitz, BP 156X, 3804-2 Grenoble France. 

(Dated: 01/03/2009) 

We have performed high resolution neutron diffraction and inelastic neutron scattering exper- 
iments in the frustrated multiferroic hexagonal compounds RMnOa (R=Ho, Yb, Sc, Y), which 
provide evidence of a strong magneto-elastic coupling in the the whole family. We can correlate the 
atomic positions, the type of magnetic structure and the nature of the spin waves whatever the R 
ion and temperature. The key parameter is the position of the Mn ions in the unit cell with respect 
to a critical threshold of 1/3, which determines the sign of the coupling between Mn triangular 
planes. 



Multiferroics have aroused a great attention for the last 
years, as the coupling between ferroelectric and magnetic 
orderings in these materials may open the route to novel 
promising electronic devices. Magnetic frustration com- 
bined with a striking magneto-elastic coupling seems to 
be at the origin of their properties, a cocktail that has a 
strong potential for novel physics These compounds 
are however rare, and far from being fully understood. 
Indeed, ferroelectricity imposes a non centrosymmetric 
space group 01, while magnetic frustration favors com- 
plex magnetic orders [3, 0] • 

The RMnOa compounds, where R is a rare earth 
ion, follow these conditions. In orthorhombic RMnOa, 
(R=Eu, Gd, Tb, Dy), where the strong GdFeOa-type dis- 
tortion lifts the orbital degeneracy, magnetic frustration 
arises from competing super-exchange interactions 
yielding incommensurate magnetic structures This 
peculiar ordering suggested a novel coupling between di- 
electric and magnetic collective modes [7]]. In hexagonal 
RMnOa with smaller R ionic radius, (R=IIo, Er, Yb, Lu, 
Y) , magnetic frustration arises from the triangular geom- 
etry yielding 120° Neel orders for the Mn moments 
0. In YMnOa and LuMnOa, where Y and Lu are non 
magnetic, an iso-structural transition was recently ob- 
served at the Neel temperature T n ^Oj , namely each ion 
"moves" inside the unit cell when Mn moments get or- 
dered. This effect provided evidence for a giant magneto- 
elastic coupling, likely connected with an increase of the 
ferroelectric polarization llj. Its origin remains unex- 
plained so far. 

To shed light on these materials, we carried out high 
resolution neutron diffraction and inelastic neutron scat- 
tering in four RMnOs (R=Y, Sc, Ho, Yb), with either 
non magnetic (Y, Sc) or magnetic (Ho, Yb) R ion, show- 
ing (Ho, Sc) or not (YjYb) a spin reorientation at Tsr 
with temperature [13, [H, [l3l- We show that the iso- 
structural transition is a systematic feature in the hexag- 
onal series. In addition, we establish a correlation be- 
tween the atomic positions, the type of magnetic struc- 
ture, and the nature of the spin waves, whatever the com- 




FIG. 1: [Colors on line]. Sketch of the hexagonal MnO 
planes with out-of-plane exchange paths Jzi (between Si 
and S3) and Jz2 (between S4, and Si as well as S4, and 
52). Black (red) arrows depict the positions of the Mn spins 
within z=0 (z=l/2) Mn planes. Si, S2 and S3 are located at 
{x,0,0), (0,a;,0), (-a;, -a;,0), while Si,SB,S6 at (x,x, 1/2), 
(1 — X, 0, 1/2) and (0, 1 — x, 1/2). Double arrows are used when 
several spin orientations can be stabilized. Inset : sketch of 
the hexagonal MnO plane with in-plane exchange paths Ji 
and J2. 



pound and its magnetic structure. We show that the key 
parameter is the position x of the Mn ions within the tri- 
angular plane with respect to a critical threshold of 1/3 
which tunes the sign of the interplane exchange interac- 
tion. We justify this result by simple energy arguments. 
Thanks to the magneto-elastic coupling, the atomic mo- 
tion helps releasing the frustration by selecting a given 
magnetic structure, depending on the x value. This pro- 
cess recalls the spin-Peierls states stabilized in several 
geometrically frustrated 2D or 3D compounds 15, lij . 
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Hexagonal RMnOs compounds consist of stacked Mn- 
O and R layers, the Mn ions forming a nearly ideal two 
dimensional triangular lattice They crystallize in 

the PG^cm space group, with two Mn-0 planes per unit 
cell. As shown on Figure [1] Mn coordinates depend on a 
unique parameter x. For x ^ 1/3, two different exchange 
path Ji and J2 can be defined between Mn moments in a 
given Mn plane (inset Fig. [1]). The triangular symmetry 
and hence the geometrical frustration within a Mn plane 
is however preserved whatever the x value. The Mn- 
Mn interactions between adjacent Mn planes are due to 
super-super exchange paths, via the apical oxygen ions 
of MnOs bipyramids. These interactions lead to a 3D 
magnetic ordering below T^r. Again, as soon as a; 7^ 1/3, 
two different paths, and thus two different interactions, 
Jzi and Jz2, can be distinguished (Fig[T]), while for x = 
1/3 all paths become equivalent. The value x = 1/3 is 
therefore a critical threshold which determines the sign 
of the effective interplane exchange Jzi — Jz2 (the same 
sign as that of a; — 1/3). As shown below, this quantity 
determines in turn the stability of the magnetic phases 
and the nature of the spin waves. 

Whatever the x value and the R magnetism, the Mn 
magnetic moments order in 120° arrangements, with 
four different possible structures (inset Fig. [2]), labeled 
from the irreducible representations (IR) F^^ i=i-4 of the 
PG^cm space group with k=0 propagation vector [ll,[lB|. 
For Fi and Mn moments are perpendicular to a and 
b axes, and their arrangement in the z = 1/2 plane are 
either antiparallel (Fi) or parallel (F4) with respect to 
z=0. The same picture holds for F2 and F3 with spins 
along a and b axes. 

High-resolution neutron powder diffraction patterns 
were collected versus temperature on the D2B and 3T2 
instruments, at ILL and LLB-Orphee reactors respec- 
tively. Powder samples HoMnOs^ScMnOa and YbMnOa 
were prepared as described in [l9| and characterized by 
x-ray diffraction. Data were analyzed using the FuUprof 
and Basireps softwares [13, [llj, allowing us to deter- 
mine the atomic positions and the magnetic structures 
versus temperature precisely. The magnetic structures 
and transition temperatures agree with previous deter- 
minations [l2| . It is worth noting that HoMnOa and 
ScMnOs undergo a second magnetic transition at Tsr 
corresponding to the reorientation of the Mn spins. 

Figure [2] shows the temperature dependence of the Mn 
position X for the three samples. As a striking feature, in 
HoMnOa, x exhibits an unprecedented large change (of 
about 3%), which occurs at Tjv, providing evidence for 
an iso-structural transition concomitant with the mag- 
netic ordering. Smaller changes occur at T/v in ScMnOa 
and YbMnOa . These results confirm those previously ob- 
tained [l3| in Y and LuMnOa, and show that this tran- 
sition is universal in the hexagonal RMnOa series. As 
noticed in Reffiol. the variations of x versus T depend of 
the rare-earth ion, namely x increases below T^r in Ho 
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FIG. 2: [Color online]. Refined Mn positions in the unit 
cell obtained from high-resolution neutron powder diffraction. 
Blue and red arrows indicate respectivly the Neel and the 
spin reorientation temperatures. Lines are guides to the eyes. 
Right insets : Mn magnetic configuration for Tsr <T < Tm- 
Left insets : Mn magnetic configuration T < Tsr). 



and YbMnOa, whereas it decreases in ScMnOa. More- 
over, we also discern important changes at the reorienta- 
tion transition. In HoMnOa, x increases with decreasing 
temperature and crosses the 1/3 threshold exactly at the 
spin reorientation temperature (Tgji). The reverse situa- 
tion holds for ScMnOa, where x decreases upon cooling, 
becoming lower than 1/3 at T^^. Finally, in YbMnOa, 
which shows no reorientation transition, x increases at 
T TV but remains lower than 1/3 in the whole tempera- 
ture range. 

All these features can be explained by considering the 
strategic position of the Mn ions. We argue that x is 
the key parameter that controls the sign of Jzi — Jz2 
and thus drives the magnetic ordering. Namely for x < 
1/3, the exchange path along Jzi is longer than along 
Jz2 and we expect Jzi — Jz2 £ 0. The reverse situation 
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TABLE I: Mn position in RMnOa compounds correlated with 
their magnetic structures defined by F irreducible representa- 
tions. YMnOs positions are taken from Ref. [l^ 



occurs for x > 1/3. A careful look at the positions and 
magnetic structures summarized in table U shows that 
X < 1/3 is associated with Fa and r4 while a; > 1/3 
corresponds to Fi and F2. This scheme also explains 
the occurrence of a spin reorientation transition when x 
crosses the 1/3 threshold. In this case, Jzi — Jz2 changes 
sign, resulting in a change from Fi towards F4 (Ho) or 
from F2 towards F3 (Sc). The crucial importance of the 
Jzi — Jz2 coupling can be further justified by calculating 
the magnetic energy £ (per unit cell) of the Mn moments 
in a mean field approximation. £ is readely obtained 
by writting the Heisenberg Hamiltonian : H — Hp + 

Hz, with Hp = J2J SiSj, and Hz — J^-^z SiSj, where 
the sums run over nearest neighbors and the subscripts 
p and z refer to in plane and out of plane interactions 
respectively. Because of the triangular arrangement, we 
have X]i=i 2 3^i ~ Si=4 5 ~ which in turn implies 



: £ = -|J5"^ + {Jzi - Jz2)iS3.Si). Depending on its 
sign, parallel or antiparallel orientations of 5*3 and S4 
are expected, giving rise to the four magnetic structures 
described above. Within this simple picture, Fi and F2 
IR minimize the energy for Jzi — Jz2 > while F3 and F4 
are favored for Jzi — Jz2 ^ 0. This is in exact agreement 
with the results summarized in Table 1. 

A straightforward way to confirm our explanation is 
to determine the value of Jzi — Jz2 by an independent 
measurement. This can be easely done by measuring the 
spin wave dispersion along the c-axis. For this purpose 
we carried out inelastic neutron scattering experiments 
on the cold triple-axis 4F spectrometer installed at LLB- 
Orphee, on large single crystals of YMn03, YbMn03 and 
HoMn03 grown by the floating zone technique. 

Figure [3I^left) shows a color map of the dynamical 
structure factor measured as a function of energy and 
wavevector in YMn03, YbMn03 and IIoMn03. These 
maps were obtained by collecting energy scans taken at 
different (1,0, Qf) wavevectors. The measurements were 
performed at 2 K in the first two cases, and at two tem- 
peratures just above and below the re-orientation tem- 
perature TsR for IIoMn03. Different features can be 
seen from these experimental data, including crystal field 
levels (Yb and Ho). A comprehensive investigation of of 
these features is however beyond the scope of this pa- 
per and we would like to focus on the low energy spin 
wave excitations labelled with arrows. We notice that 



in YMn03 and YbMn03 , this particular branch displays 
upwards (Yb) or downwards (Y) dispersions, revealing 
opposite couplings along c. Similarly in HoMn03, when 
crossing the reorientation "TsR, the curvature changes 
from upwards to downwards, indicating a reversal of the 
magnetic interaction along the c-axis. 

To get a quantitative information about these cou- 
plings, we performed a spin wave analysis of the Heisen- 
berg Hamiltonian H = Hp + Hz defined above, taking 
into account additional planar and uniaxial anisotropy 
26|. Figure [ll^right) shows the dynamical 
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terms 

structure factors calculated on the basis of this model as 
a function of energy transfer w and wavevector (1, 0, Qi). 
In this approach, six spin wave modes are expected. 
Along c, four of them are almost degenerate. They ex- 
hibit a large gap due to the planar anisotropy, as well as 
a very weak dispersion. On Figure [3] (right), these modes 
correspond to the flat branch sitting at 5 meV (Ho,Y) or 
6 meV (Yb). The two remaining modes correspond to 
the Goldstonc modes of the magnetic structure. In g = 
limit, they can be seen as global rotations of the 120° 
pattern inside the basal planes. They couple via the Jz 
interactions, resulting in either in phase or out of phase 
rotations. These two modes are pointed out by arrows 
in Figure [3] (right) . Actually, one of them has a van- 
ishing intensity but can still be observed around (101). 
Assuming antiferromagnetic couplings Jzi and Jz2, we 
can determine from the data Jzi — Jz2 — 0.0050(5) meV 
in YMn03 and Jzi - Jz2 = -0.012(2) meV in YbMn03. 
For HoMn03, the neutron data are well modelled with 
Jzi - Jz2 = -0.0038(5) meV at T=45K > Tsr and with 
Jzi - Jz2 = 0.0018(5) meV at T=27K < Ts^. We note 
that the corresponding signs of Jzi — Jz2 deduced from 
these measurements are fully consistent with the diffrac- 
tion analysis. 

Our analysis emphasizes the double origin of the mag- 
netic frustration. The Mn triangular planes are geomet- 
rically frustrated for antiferromagnetic interactions. In 
addition, adjacent Mn planes are coupled along the c 
axis by self-competing interactions. The Mn shift with 
respect to the 1/3 position does not suppress the rota- 
tional invariance in the Mn plane, but clearly lifts the 
interplane frustration, allowing 3D ordering. Depending 
on X, either (Fi,F2) or (F3,F4) structures are stabilized. 
The remaining degree of freedom in the system is the 
global rotation of a 120° Neel order around the c-axis. 
Namely, to stabilize Fi or F2 (F3 or F4), the Mn mag- 
netic moments must couple to the (a, 5) crystal axes. 

On the one hand, residual anisotropic interactions 
among Mn ions can in principle select a given orienta- 
tion in a Mn plane. As revealed by the exceptionally 
small uniaxial gap observed in YMn03 2^, 2^ , these in- 
teractions are weak. They are not sufficient to break the 
triangular symmetry, so that two-dimensional spin liquid 
ffuctuations remain |22| . On the other hand, one could 
argue that R-Mn interactions (when R is magnetic) play 
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FIG. 3: (color online.) Maps of the dynamical structure fac- 
tor in Ho, Yb and YMnOs. The neutron intensity is plotted 
versus the energy transfer lu and wave vector (1, 0, Qi). Left 
panel : experimental results from inelastic neutron scatter- 
ing. Black lines highlight the low energy spin wave modes. 
Additional Q-independent lines in Yb and HoMnOs maps are 
due to crystal field excitations. Right panel : spin wave dis- 
persion curves calculated using an Heisenberg Hamiltonian. 
Black arrows indicate the two Goldstone modes (see text). 



a significant role. Indeed, the R moments on the 46 site 
order at Tjv in the Mn molecular field, while their orien- 
tations are clearly coupled to the Mn ones, through an 
energy term of anisotropic nature Our spin wave 

measurements (Fig[3|) show that the uniaxial anisotropy 
gap in YbMnOa (3 meV) varies with temperature like the 
Yb moment, providing evidence for such R-Mn coupling. 
Nevertheless, we shall not conclude that the Mn orienta- 
tion in the (a, b) plane is completely determined by this 
interaction. For instance, spin reorientation transitions 
of the Mn sublattice may occur (Sc, Ho) or not (Y, Yb) 
whatever the R magnetism. 

In conclusion, our elastic and inelastic neutron scatter- 



ing experiments clearly show the crucial role of the Mn 
position in determining both the magnetic structure and 
the spin waves modes. The onset of Mn magnetic order- 
ings at Tjv or Tsr correlates with the Mn position. The 
magnetic orders and spin excitations in the whole series 
result from a subtle interplay of magneto-elastic coupling, 
frustrated intra- and inter-plane Mn-Mn interactions and 
R-Mn interactions. The strong importance of inter-plane 
interactions strongly suggests that the shift of the Mn 
position, which releases the frustration along the c axis, 
is a key ingredient at the origin of the multiferroicity in 
the hexagonal RMnOa family. We thank E. Suard for 
her help in the D2B experiment and V. Simonet and F. 
Damay-Rowe for useful discussions. 
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